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ABSTRACT 

We present a multiwavelength study of the infrared dark cloud MSXDC G034. 43+00. 24. Dust emis- 
sion, traced by millimeter/submillimeter images obtained with the IRAM, JCMT, and CSO telescopes, 
reveals three compact cores within this infrared dark cloud with masses of 170-800 Mq and sizes < 
0.5 pc. Spitzer 3.6-8.0 /im images show slightly extended emission toward these cores, with a spec- 
tral enhancement at 4.5 /zm that probably arises from shocked H2. In addition, the broad line widths 
(AV ~ lOkms- 1 ) of HCN (4-3) and CS (3-2), and the detection of SiO (2-1), observed with the 
JCMT and IRAM telescopes, also indicate active star formation. Spitzer 24 /jm images reveal that 
each of these cores contains a bright, unresolved continuum source; these sources are most likely em- 
bedded protostars. Their millimeter to mid-IR continuum spectral energy distributions reveal very high 
luminosities, 9000-32,000 Lq. Because such large luminosities cannot arise from low-mass protostars, 
MSXDC G034.43+00.24 is actively forming massive (~10 M Q ) stars. 

Subject headings: infrared: stars-ISM: molecules-stars: formation-submillimctcr 



1. INTRODUCTION 

Infrared dark clouds (IRDCs) are a distinct class of in- 
terstellar gas cloud discovered in the Infrared Space Obser- 
vatory and Midcourse Space Experiment (MSX) Galactic 
plane mid-IR surveys (Perault et al. 1996; Carey et al. 
1998; Hennebelle et al. 2001). IRDCs are identified as dark 
extinction features seen in silhouette against the bright 
Galactic background at mid-IR wavelengths. Preliminary 
studies show that they have high densities (> 10 5 cm~ 3 ), 
high column densities (~ 10 23 — 10 25 cm -2 ), and low tem- 
peratures (< 25 K; Egan et al. 1998; Carey et al. 1998, 
2000). 

It is still unclear what role IRDCs play in star forma- 
tion. Egan et al. (1998) suggested that IRDCs are quies- 
cent, isolated molecular clumps, not yet actively forming 
stars. Indeed, if a massive star formed inside an IRDC, it 
would quickly heat the surrounding natal dust, and thus 
the IRDC would not remain cold and dark. The detection 
of compact millimeter and submillimeter (Lis & Carlstrom 
1994; Carey et al. 2000; Redman et al. 2003) continuum 
sources toward some IRDCs, however, suggests that they 
may harbor prcstcllar cores. If so, then IRDCs may well 
play a key role in the very earliest stages of star formation. 
It is important to determine whether these compact cores 
are sites of current or future star formation. 

In this Letter, we use two standard methods to iden- 
tify sites of active star formation. The first relies on the 
interaction between a young stellar object (YSO) and its 
parental molecular cloud. Because YSOs undergo bipo- 
lar outflows, they interact strongly with the surrounding 



gas and thus can be found by detecting large molecular 
line widths or shocked gas. One can search for shocks 
by detecting spectral lines that require shocks for exci- 
tation (e.g., H2) or chemical species that require shocks 
for their formation (e.g., SiO). The second method relies 
on the direct detection of the protostar. Because proto- 
stars are deeply embedded in dusty clouds, they emit most 
strongly in the millimeter to IR. Active star-forming cores 
are thus associated with luminous, highly reddened, com- 
pact sources. 

We have used both of these methods to discover 
three distinct sites of active star formation in the IRDC 
MSXDC G034.43+00.24 (hereafter G34.4+0.2). We iden- 
tified this IRDC as an extinction feature in the MSX 8 /im 
Galactic plane survey (Simon et al. 2005). Based on its 
13 CO emission from the BU-FCRAO Galactic Ring Sur- 
vey (Simon et al. 2001), we find that G34.4+0.2 is located 
at a kinematic distance of 3.7 kpc (because it is an IRDC, 
we assume it lies at the near kinematic distance). 

Superposed against this molecular cloud is the compact 
H 11 region IRAS 18507+0121 (Bronfman et al. 1996; Moli- 
nari et al. 1996). Shepherd et al. (2004) noted a millimeter 
continuum source (our core MM1) ~40" to the north of 
this Hn region. From its weak 6 cm emission, they sug- 
gest that this millimeter core contains a deeply embedded 
B2 protostar. Although Shepherd et al. (2004) did not 
note that this millimeter core is in fact part of an IRDC, 
we find that it is indeed embedded within G34.4+0.2. In 
addition, Garay et al. (2004) detected 1.2 mm continuum 
emission associated with this cloud in their survey of mas- 
sive star-forming regions. They note that one core (our 
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MM3) appears to be a cold (<17 K), massive, dense, and 
quiescent core. 

We find strong evidence to support the idea that these 
two cores, and one other within G34.4+0.2, are indeed 
sites of active star formation and will eventually become 
massive stars (M ~ 10 M©). The discovery of active mas- 
sive star formation in this IRDC suggests that IRDCs may 
be the earliest stage in the formation of high-mass stars. 

2. OBSERVATIONS 

We present both continuum and spectral line data 
at many wavelengths (see Table 1 for details). The 
continuum data cover the wavelength range from 1.2 
mm to 24 fim. The millimctcr/submillimctcr continuum 
data were taken at the Institute de Radioastronomie 
Millimetrique (IRAM), James Clerk Maxwell Telescope 
(JCMT), and Caltech Submillimeter Observatory (CSO) 
using either an "on-the-fly" or "scan-mapping" mode. 
Standard reduction methods within the software pack- 
ages MOPSIC, CRUSH, and SURF were used to correct 
for atmospheric opacity and to remove atmospheric fluc- 
tuations. The millimetcr/submillimetcr data were flux- 
calibrated using either G34.3 or Uranus. The IR con- 
tinuum data were obtained using the Spitzer Space Tele- 
scope. The 3.5-8.0/zm images were obtained as part of 
the Galactic Legacy Infrared Mid-Plane Survey Extraordi- 
naire (GLIMPSE; Benjamin et al. 2003) using the Infrared 
Array Camera (IRAC). The 24 /im images were obtained 
using the photometry raster mapping mode of the Multi- 
band Imaging Photometer for Spitzer (MIPS). The IRAC 
data were processed by the GLIMPSE team, and the MIPS 
data were processed by the Spitzer Science Center data 
processing pipelines. 

The spectroscopic data consist of millime- 
ter/submillimctcr molecular line observations obtained 
with IRAM and JCMT. Molecular line spectra toward two 
of the millimeter cores, MM1 and MM3, within G34.4+0.2 
were obtained in SiO (2-1), CS (3-2), HCN (4-3), and 
C 18 (3-2). Chopper wheel calibration, pointing, and 
focus checks were performed regularly. The data were 
reduced using standard methods in the software packages 
GILDAS and SPECX. 

3. RESULTS AND DISCUSSION 

3.1. Identifying Candidate Protostellar Cores 

To locate potential sites of star formation, we have im- 
aged G34.4+0.2 in the millimeter/submillimcter contin- 
uum. Because these wavelengths probe cold dust emis- 
sion, millimcter/submillimeter continuum emission can lo- 
cate dense compact cores, potential sites of star formation. 
Our millimetcr/submillimetcr images of G34.4+0.2 reveal 
extended, diffuse continuum dust emission that matches 
the morphology of the mid-IR extinction extremely well. 
In addition to this extended emission, we also find four 
compact cores. Figure 1 shows the 1.2 millimeter con- 
tinuum emission overlaid on a Spitzer three-color image 
(3.6 /im in blue, 4.5 /im in green, and 8.0 /im in red). One 
of these cores (MM2) is clearly associated with the Hn 
region IRAS 18507+0121. The remaining three dust cores 
(MM1, MM3, and MM4) are unresolved (< 0.5 pc) and 
have 1.2 mm peak fluxes that range from 0.3 to 2.5 Jy 
beam -1 . Table 2 lists their coordinates and peak fluxes. 



Because the 1.2 mm continuum emission is optically 
thin, we can reliably estimate masses via the expression 
M = F„D 2 /[k 1/ B 1/ (T)}, where F„ is the observed flux den- 
sity, D is the distance, k v is the dust opacity per gram of 
dust, and B„(T) is the Planck function at the dust tem- 
perature. We adopt for «i.2mm a value of 1.0 cm 2 g _1 
(Ossenkopf & Henning 1994). Integrating the 1.2 mm con- 
tinuum emission above a 3 a level, and assuming a dust 
temperature of = 30 K and a gas-to-dust mass ratio of 
100, we find that the total mass of the IRDC is -7500 M Q . 
Gaussian fits to the three cores reveal masses in the range 
170-800 M© (see Table 2). The core masses and sizes 
are similar to massive protostars seen toward other star- 
forming regions (Brand et al. 2001). 

3.2. Evidence of Outflows and Shocked Gas 

Evidence of kinematic interactions between embedded 
YSOs and the molecular gas in these cores can be found 
in the Spitzer/IRAC data. The IRAC three-color 3.5, 
4.6, and 8.0 [im image of G34.4+0.2 reveals slightly ex- 
tended emission toward all three millimeter cores, with a 
relative enhancement at 4.5 /mi (displayed as green). Al- 
though highly extincted stars can appear green in these 
three color images due to the flat extinction curve between 
4.5 and 8.0/im (Indebetouw et al. 2005), bright, extended 
4.5 /um emission has been found to trace shocked gas in 
star-forming regions (Marston et al. 2004; Noriega-Crespo 
et al. 2004). This enhanced 4.5 /um emission is thought 
to arise from shock-excited spectral lines, probably the H2 
5(9) line at 4.694 /im. Because these cores show evidence 
of shocked gas, they are probably forming stars. 

Further evidence comes from our IRAM and JCMT 
molecular line spectra toward two of these cores (Fig. 2). 
The high-density tracing molecular lines, CS (3-2) and 
HCN (4-3), show broad line widths (AV -10 kms" 1 ). 
Such broad line widths are often interpreted as tracing pro- 
tostellar outflows (Moriarty-Schieven et al. 1995). Because 
SiO is thought to form as the result of the shock destruc- 
tion of dust grains, the SiO (2-1) line indicates strong 
shocks and is therefore a reliable indicator of star for- 
mation (Bachiller et al. 1991; Martin-Pintado et al. 1992; 
Zhang et al. 1995; Gibb et al. 2004). The strong SiO line, 
combined with broad lines of the molecular high-density 
tracers, provides further evidence that active star forma- 
tion is occurring within these cores. 

3.3. Embedded Protostars 

A more direct indicator of active star formation is the 
detection of embedded protostars via their millimeter to 
IR continuum emission. We have detected embedded pro- 
tostars toward the three cores using new Spitzer 24 /im 
data. The 24 /im image of G34.4+0.2 shows that, while 
the IRDC itself remains dark, the active cores each con- 
tain unresolved (<6"), bright 24 /im continuum sources 
(Fig. 1). 

We can determine the luminosity of the embedded pro- 
tostars from their millimeter to IR broadband spectral en- 
ergy distributions (SEDs). Using the peak fluxes from all 
of our continuum data in the mid-IR (24 /im), submil- 
limeter (350, 450, and 850 /im), and millimeter (1.2 mm), 
we can construct broadband SEDs and fit them to stan- 
dard graybody curves (Fig. 3). For simplicity, we assume 
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that the cores are isothermal. In our fits, the emissiv- 
ity index (/?), optical depth at 250 /xm (r 2 5o), and dust 
temperature (To) are free parameters. For a source size 
of 15", these fits yield reasonable values for these pa- 
rameters: /3~1.8, r 25 o~0.25, and Td~33K (see Table 2). 
The SEDs also provide good estimates for the bolometric 
luminosity, which is simply the integral under the SED 
fit. Because these sources are slightly resolved, fits to 
the SEDs using the peak fluxes will yield lower limits to 
the actual luminosity. We find that the three protostcllar 
cores in G34.4+0.2 have bolometric luminosities (Lf, i) of 
9000-32,000 L Q . 

3.4. Large Luminosities: High-Mass Protostars 

Theoretical studies of protostcllar evolution suggest 
that low-mass stars (M<2 M Q ) never achieve luminosi- 
ties > 100 Lq in their prc-main-sequcnce evolution (Ibcn 
1965). High-mass protostars (M>8 Mq), however, evolve 
at essentially constant luminosity on their path from pro- 
tostar to main-sequence star (Iben 1965; Palla & Stahlcr 
1990). Thus, the protostcllar luminosity is a good indica- 
tor of the eventual main-sequence luminosity and, hence, 
stellar mass. Protostellar luminosities of — 10,000 Lq, 
therefore, correspond to early B stars with masses of 
-10 M Q (Panagia 1973). Because G34.4+0.2 contains 
three sites of active star formation with Lt, i —10,000 L & 
or more, we conclude that each of these sites contains 
a massive protostar. Specifically, wc find corresponding 
spectral types of 09.5 for MM1 and B0. 5 for both MM3 
and MM4 (Panagia 1973). 

3.5. The Role of IRDCs in Star Formation 

The discovery of active massive star formation in 
G34.4+0.2 shows that the earliest stages of massive star 
formation may well take place in IRDCs. Because low- 
mass stars form within Bok globules, high-mass stars prob- 
ably form within a higher mass equivalent. Wc speculate 
that IRDC cores are indeed these high-mass counterparts 
to Bok globules. If IRDCs host the earliest stages of high- 
mass star formation, and then at some point they will 
spawn OB stars and grow warm. Indeed, the masses and 
sizes of IRDCs (M~ 1000 Mq, sizes approximately a few 
parsecs) are similar to those of warm molecular clumps 
associated with massive star-forming regions (Garay & 
Lizano 1999; Lada & Lada 2003). Moreover, such warm 
clumps often produce multiple massive stars. The fact that 
G34.4+0.2 is forming three massive stars and has already 
formed a B0. 5 star associated with the Hn region IRAS 



18507+0121 (Shepherd et al. 2004) suggests that IRDCs 
can also form several massive stars almost simultaneously. 
To test this hypothesis, we are currently searching for more 
examples of active star formation within a larger sample 
of IRDCs. 

4. CONCLUSIONS 

In order to search for regions of active star forma- 
tion in IRDCs, wc have performed a multiwavelcngth 
study of the IRDC G34.4+0.2. We identified three po- 
tential star-forming cores by finding compact millime- 
ter/submillimctcr continuum sources. These cores show 
enhanced 4.5 /im emission, which is thought to trace 
shocked gas in star-forming regions. Broad molecular line 
widths of CS and HCN provide additional evidence of star 
formation toward two cores. Each of these cores also has 
strong, broad lines of SiO, a species that preferentially 
forms in shocks. 

We have also directly detected bright, highly reddened, 
compact sources, presumably protostars, in all three cores. 
The broadband SEDs of these cores, from millimeter to 
mid-IR wavelengths, combined with the known kinematic 
distance, provide a good estimate of their bolometric lumi- 
nosities, 9000-32,000 L Q . Both theoretical and empirical 
studies suggest that such large luminosities can only arise 
from high-mass protostars. Hence, G34.4+0.2 contains 
three active sites of star formation, which will eventually 
produce massive stars with M — 10 M Q . This study sug- 
gests that IRDCs play an important role in the early stages 
of massive star formation. 
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Table 1 
Summary of observations 



Telescope 


Instrument 


Wavelength/ 


Date 


Angular 


lcr rms 






Spectral line 




Resolution (") 


noise 


IRAM 30 m 


MAMBO-II 


1.2 mm 


2004 Feb 


11 


10 mjy 


JCMT 15 m 


SCUBA 


850, 450 ^m 


2004 Sept 


15, 8 


100, 300 mjy 


CSO 10 m 


SHARC-2 


350 /im 


2005 Apr 


8 


500 mjy 


Spitzer 


MIPS 


24 /mr 


2004 Oct 


6 


70 ^Jy 




IRAC 


3.6, 4.5, 8.0/im 


2003-2004 


2 


23, 27, 78 ^Jy 


IRAM 30 m 


C150 Rx 


CS (3-2) 


2004 Nov 


26 a 


0.13 K 




A100 Rx 


SiO (2-1) 


2004 Nov 


28 


0.06 K 


JCMT 15 m 


B3 Rx 


HCN (4-3), C 18 (3-2) 


2004 Apr 


14 


0.07, 0.13 K 



"Mapped and then smoothed to this resolution 
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Table 2 

Parameters for the 1.2 mm cores within G34.4+0.2. 



Source a 


Coordinates 






Peak Flux 










SED fit parameters 




a S 


1.2mm 


850^m 


450^m 


350/im 


24^111 


Mi.2mm b 





T250 Tj) Ltol 




(J2000) (J2000) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(M ) 




(K) (L Q ) 


MM1 


18 53 18.0 01 25 24 


2.5 


6.6 


70 


140 


>0.09 


800 


U 


, 0.53 34 32,000 


MM3 


18 53 20.4 01 28 23 


0.3 


1.2 


10 


20 


0.03 


200 


u 


, 0.08 32 9,000 


MM4 


18 53 19.0 01 24 08 


0.6 


1.8 


15 


30 


0.04 


170 


li 


, 0.13 32 12,000 



a The cores are named in order of peak 1.2 mm flux (determined from gaussian fits). IRAS 18507+0121 corresponds to our core MM2. 
^Calculated from a gaussian fit to the 1.2 mm emission associated with the core (excluding the underlying emission from the IRDC). 




Right Ascension (J2000) Right Ascension (J2000) 

Fig. 1. — G34.4+0.2. Left: Spitzer/IRAC three-color image (3.6 ^m in blue, 4.5 fim in green and 8.0 /an in red) overlaid with 
IRAM /MAMBO-II 1.2 mm continuum emission (contour levels 60, 90, 120, 240, 480, 1200, and 2200 mjy beam -1 ). Right: Spitzer/MIPS 
24 fira image with contours of the IRAM/MAMBO-II 1.2 mm continuum emission (logarithmic color scale, 30 Mjy sr _1 [black] to 200 Mjy 
sr _1 [white]). Labelled on this figure are the four millimeter cores. 



6 



Rathborne et al. 




40 50 60 70 80 40 50 60 70 80 
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Fig. 2. — IRAM and JCMT molecular line emission toward two of the millimeter cores. The solid vertical line marks the central velocity 
of the core [as traced by the optically thin C ls O (3-2) emission]. Note the broad line widths in the HCN (4—3) and CS (3—2) spectra. 
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Fig. 3. — Broadband continuum SED for the three cores. Included on this plot are peak fluxes at 24 /im, 350 fim, 450 fim, 850 (im, and 
1.2 mm. The curves are graybody fits to the data, which yield values of /3~1.8, T25o~0.25, and Trj~33K (15" source diameter; see also 
Table 2). The derived bolometric luminosities are labeled for each core. MM1 saturates the MIPS array; hence, the quoted 24 (im flux is a 
lower limit. 



